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Abstract 

In the framework of the scale relativity theory, the chaotic be- 
havior in time only of a number of macroscopic systems corresponds 
to motion in a space with geodesies of fractal dimension 2 and leads 
to its representation by a Schrodinger-like equation acting in the 
macroscopic domain. The fluid interpretation of such a Schrodinger 
equation yields Euler and Navier-Stokes equations. We therefore 
choose to extend this formalism to study the properties of a system 
exhibiting a chaotic behavior both in space and time which amounts 
to consider them as issued from the geodesic features of a mathe- 
matical object exhibiting all the properties of a fractal 'space-time'. 
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Starting with the simplest Klein-Gordon-hke form that can be given 
to the geodesic equation in this case, we obtain a motion equation 
for a 'three fluid' velocity field and three continuity equations, to- 
gether with parametric expressions for the three velocity components 
which allow us to derive relations between their non- vanishing curls. 
At the non relativistic limit and owing to the physical properties 
exhibited by this solution, we suggest that it could represent some 
kind of three-dimensional chaotic behavior in a classical fluid, tenta- 
tively turbulent if particular conditions are fulfilled. The appearance 
of a transition parameter I? in the equations allows us to consider 
different ways of testing experimentally our proposal. 

KEY WORDS: Fluid mechanics; Chaotic behavior; Scale relativ- 

ity. 

1 Introduction 

The Madelung transformation of the Schrodinger equation has been 
long known as providing a tentative fluid-like interpretation of non 
relativistic quantum mechanics [TJ [2]. When writing the complex 
wave function in polar form, t/j — \/Pe'*, then separating the 
real and imaginary parts of the Schrodinger equation, one obtains a 
Hamilton- Jacobi and a continuity equations where P and V = \79 
can be interpreted as, respectively, the density and the velocity field 
of an irrotational flow in a compressible fluid. The gradient of the 
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Hamilton- Jacobi equation is an Euler-type equation exhibiting a new 
term, the so-caUed 'quantum' potential, which depends only on the 
density and has been, in the hydrodynamical picture, regarded as a 
pressure potential or a mechanical stress tensor [SlU]. 

In scale relativity, the effect on motion of the simplest scale laws, 
namely self-similar laws with fractal dimension 2, allows one to re- 
cover the various equations of quantum mechanics [5] [6l [71 18] . How- 
ever, since the mathematical construction of the quantum equations 
thus obtained does not depend on the value of the Planck constant 
h (see Sec. |3]), its results can be generalized to the classical realm 
provided some particular conditions are fulfilled by the system under 
consideration 5, 7 . This has already been done for the Schrodinger 
equation with applications to the chaotic behavior in time of a num- 
ber of macroscopic systems (planetary systems [5l[9l[T0], astrophys- 
ical structures [111 I12j . etc.) which have been successfully corrobo- 
rated by observations. 

Now, it must be remembered that it was shown that non-differentiability, 
which is at the root of the theory, implies a fractal space(-time) 
which in turn implies an infinite number of geodesies defined at any 
of its points, hence a fiuid-like interpretation of such geodesic bun- 
dles O [H 1131 114] . In the case where the fractality of space only is 
involved, the parameter along the geodesies is the time t and the 
differentiable elements representing the resolutions which define the 
scale state of the system under consideration is dt. A fluid interpre- 
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tation of the macroscopic Schrodinger-type equation can therefore 
be developed which allows one to obtain classical fluid motion equa- 
tions, i. e., Euler and Navier-Stokes equations, with a 'quantum- like' 
potential [15j . 

To strengthen the appropriateness of the use of scale relativity 
methods to deal with chaotic motion, we stress that the quantum- 
type equations, from which the fluid motion equations are derived, 
are obtained in this framework when three conditions are fulfilled. 
They are consequences of the giving up of the differentiability as- 
sumption and read, in the simplest case when only the dt ^ — dt 
symmetry is broken [5]: 

1. The geodesic number is infinite as a consequence of the frac- 
tality of space. Thus, a fluid-like description obtains in which the 
velocity v{t) is replaced by a velocity field v[x{t),t]. 

2. The geodesies are fractal curves with fractal dimension 2 which 
is the dimension of a Markovian process that we take as a represen- 
tation of the assumed chaotic behavior of our system. 

3. Due to non-differentiability, the invariance under the differ- 
ential element symmetry dt ^ — dt is broken, in the sense that the 
standard definition of a derivative where dt is assumed to vanish 
does no more apply. One therefore replaces the coordinates x{t) by 
fractal functions x{t, dt) and the velocity becomes two- valued as 



v+ [x{t, dt),t, dt] 



x{t + dt, dt) 



x{t, dt) 



(1) 



dt 
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V- [x(t, di),t, dt\ 



x{t, dt) — x{t — dt, dt) 
Jt 



(2) 



This yields two fractal velocity fields 



V+ = v+[x{t),t\ + w+[x{t,dt),t,dt\, 



(3) 



V- = v-[x{t),t\ + w-[x[t,dt),t,dt\, 



(4) 



where v+, respectively u_, is a 'classical part', differentiable and 
independent of the resolution dt, and ■w+, respectively W-, is a 'frac- 
tal part', explicitly dependent on dt and divergent when dt goes to 
zero. The 'classical part' is shown to dominate over the 'fractal part' 
above (or below) some transition scale characterized by a parameter 
D [5l E], hence two different behaviors of the system: non-chaotic 
beyond (below) this scale and chaotic below (beyond). 

When one wants to go farther on, one must consider that the 
parameter along the geodesies is the proper time s and thus, as 
a first step, consider the breaking of the space-time differential el- 
ement ds, which, when considered alone, gives the complex Klein- 
Gordon equation [7] . This amounts to analyze only the effect of non- 
differentiability on the total derivative d/ds. Ifowever, the velocity 
fields of the 'geodesic' bundles are also functions of the coordinates, 
so that one is led to consider also the effect of the partial derivatives 
d/dx^ in the decomposition d/ds = d/ds -f [dx^ / ds)d / dx^ which 
yields the new symmetry breaking dx^ ^ ~ dx'^ [6]. In this case, 
the corresponding algebra is that of the quaternions. 

Note that, since we are now in a full space-time representation. 
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the time t, which was a mere parameter on the geodesies in the pre- 
vious case, becomes here a coordinate of a fractal four-dimensional 
space-time. Thus, while projecting the fractal geodesic curves on 
such a coordinate, we see that time becomes also a fractal function, 
as it is the case of the three spatial coordinates. 

Therefore, we propose here to examine whether the chaotic be- 
havior both in time and space of the macroscopic system we con- 
sider can be constructed as emerging from the breaking, at the level 
of an underlying mathematical object which have all the properties 
of a fractal 'space-time', of both the differential element symmetries 
ds ^ — ds and dx^ ^ — dx^ . The development of such a formalism, 
which yields a quaternionic Klein-Gordon equation, is followed by a 
transition to the non-relativistic limit, since our aim is to model non- 
relativistic flows. As a first step, we limit ourselves to the simplest 
motion equation which can be written in this framework and which 
corresponds to free (inertial) motion in the fractal 'space-time'. We 
leave the study of more complicated behavior to future work. How- 
ever, we are thus able to show that new and interesting properties 
emerge as well as a new way of characterizing the possible transition 
to some kind of 'turbulence'. 

A long-known property of the quaternionic formalism is that the 
Dirac equation for a free particle can be obtained as a mere square 
root of the Klein-Gordon operator [161 117) . This property has been 
used to derive the Dirac equation as a geodesic equation in a non- 
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difFerentiable therefore fractal space-time in the case when the ve- 
locity field of the geodesic fiuid, and therefore the wave function, is 
represented by bi-quaternions (i. e., complex quaternions) [51 113], 
The bi-quaternionic wave function has been shown to exhibit all 
the properties of Dirac bi-spinors. Moreover, when one takes the 
non-relativistic limit of this bi-quaternionic equation, one recovers 
the Pauli equation and all the well known properties of the non- 
relativistic Pauli spinor [8]. We therefore suspected, in a first stage 
of our study, that relativistic fluids whose motion equations could be 
obtained by such a process might have richer properties than those 
of the simple fluids obtained in the complex wave function case. 

However, when trying to use this method in order to obtain some 
chaotic behavior for a classical fluid as the non-relativistic limit of a 
relativistic bi-quaternionic one, the vorticities of the velocity fields 
disappear and we are left with a set of gradients which are not well 
designed to represent a whirling fluid. 

This comes from the way algebra doublings appear in this con- 
text (see Sec. [21. Actually, to each symmetry breaking, ds ^ — ds, 
dx^ ^ —dx^, then x^ 'r-* ~ corresponds a doubling of the alge- 
bra, which allows one to use, first, the complex numbers, then the 
quaternions and finally the bi-quaternions, and their related algebras 
[6]. However, the breaking of the symmetry x^ —x^ is a, specific 
property of the microscopic quantum domain which is not observed 
in classical physics. Therefore, we do not need to consider here the 
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breaking of this symmetry and we are left with only two symmetry 
breakings and therefore with the quaternionic algebra. 

A first interesting attempt was made in this context by Love and 
Boghosian flSl. However, these authors consider a Madelung trans- 
formation of the quaternionic Schrodinger equation without giving 
any physical justification for this choice. Only a mathematical con- 
sideration is put forward to explain it. Note moreover that, in the 
scale relativity framework, the Schrodinger equation can be obtained 
from only one symmetry breaking, namely that of the differential 
time element dt. Therefore its corresponding algebra is that of the 
complex numbers, not of the quaternions. 

The content of this paper is as follows. In Sec. [51 we give a re- 
minder of the method which allows one to recover the Klein- Gordon 
equation with a complex wave function and the standard Euler equa- 
tion in the context of scale relativity. This is provided in order to 
give the reader a simple example of how the scale relativity methods 
can be used to obtain classical motion equations of fluid mechan- 
ics from quantum-like mechanical equations assumed to act in the 
macroscopic domain. In Sec.O we jump to the quaternionic case and 
derive a set of motion equations for the three 4-vector components 
of a 'relativistic' non-Abelian fluid velocity fleld, together with para- 
metric expressions for these components which exhibit their main 
physical properties. In Sec. |4l we take the non-relativistic limit of 
the solution obtained in Sec. |3] which gives motion equations for the 
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fluid velocity field which appears under the form of a set of three 
vector fields. Three conservation equations and expressions for the 
non vanishing curls of these three vector fields are also obtained at 
this limit. Owing to the way we have constructed this solution from 
a scale relativistic treatment of chaos, we propose that, provided par- 
ticular conditions are fulfilled, some kind of chaotic motion, possibly 
turbulent, might be represented by such a 'three fluid' . A discussion 
and our conclusion are given in Sec. \5\ 



2 Complex Klein-Gordon equation and 
Euler equation 

Following Madelung, some authors have applied the fluid representa- 
tion to the complex Klein-Gordon field of relativistic quantum me- 
chanics O 131 [19]. They also obtain an Euler- type equation for a 
relativistic fiuid and a conservation equation for a relativistic current- 
type quantity j'^ . However, one can derive directly the fluid motion 
equations from geodesic equations written with the scale relativity 
covariant derivative without any need to integrate them previously 
under the form of a Klein-Gordon equation (19) . This is the method 
that we apply here. 

The example of a scale relativity method implementation given 
in this Section is provided as a summarized reminder. Since a fully 
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detailed development and discussion of all the aspects of the formal- 
ism is beyond the scope of this paper, we refer the interested reader 
to the original articles. 

Scale relativity generalizes to scales the relativity principle only 
applied to position, orientation and motion in other standard rela- 
tivity theories. 

Here, space-time is assumed to be a continuum which is not only 
curved, as in general relativity, but can also be non-differentiable 
and therefore fractal [5]. The infinite number of fractal geodesies at 
each point of such a fractal space-time allows a fluid representation 
of the theory (see [M] for a discussion). The effect on motion of the 
simplest scale laws, namely self-similar laws with fractal dimension 2, 
allows one to recover the various equations of quantum mechanics [6] 
[71 . They are derived as geodesic equations each corresponding to 
a given level of discrete symmetry breaking, namely of the reflection 
invariances of differential elements, which are implicit symmetries of 
classical mechanics. 

The Schrodinger equation is obtained as the integral of a geodesic 
equation acting in a fractal space and issued from the breaking of the 
reflection invariance of the time differential element dt ^ — dt 
The Klein-Gordon equation proceeds from a geodesic equation acting 
in the fractal space-time and from the breaking of the relativistic line 
element symmetry ds ^ — ds [7] [20]. In both cases the wave function 
is complex as a consequence of the two-valuedness issued from these 
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symmetry breakings. 

Then, a new two-valuedness is introduced, in the motion-relativistic 
case, from the breaking of the discrete symmetry of the space-time 
coordinate differential elements dx^ ^ — dx^ . It is naturally de- 
scribed in terms of quaternions (spinors). Adding parity violation, 
i. e., the breaking of the discrete symmetry x^ ^ — , which is 
a property of standard quantum theory, one is led to introduce bi- 
quaternions (bi-spinors) . The Dirac equation is then derived as the 
bi-quaternionic square root of a bi-quaternionic Klein-Gordon equa- 
tion issued from a geodesic equation acting in the fractal space-time 
[H [13]. The Pauli equation is recovered as the quaternionic non- 
relativistic limit of this Dirac equation [S] . 

The fluid interpretation of the Schrodinger equation in the frame- 
work of scale relativity was developed in [15II14|[T9] . A generalization 
to the relativistic complex Klein-Gordon case was given in [19]. We 
recall briefly below the main steps of the reasoning since it will be 
adapted to the quaternionic case in Sec. [S] Instead of performing 
a Madelung-type transformation of the Klein-Gordon equation, the 
idea is to take advantage of the formulation of this equation as a 
geodesic equation in the fractal space-time to obtain the fluid mo- 
tion equation directly from the covariant equations of the dynamics. 

We have seen in Sec. [T] that, in the scale relativity approach, 
the velocity field on the geodesic bundle of the non-differentiable 
space(-time) is decomposed in terms of a differentiable part, V, and 
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of a fractal, divergent, non-differentiable part, W (of zero mean). 
At this stage of the construction, both are complex due to the fun- 
damental two-valuedness issued from the breaking of the symmetry 
ds ^ — ds. A complex covariant total derivative can be built from 
the different iable part of the velocity field V which, when applying 
to the relativistic case, is [B] HD] 



where the constant A gives the amplitude of the fractal fluctuations 
d^, as A = {d^^)/ds (in the quantum mechanics interpretation, A is 
the Compton length of the 'particle', A = h/mc). Now one assumes 
that the system can be characterized by an action S which is complex 
such as the four- velocity. Following the same reasoning as in classical 
mechanics, one writes 



The stationary action principle applied on this action yields a free 
equation of motion which takes the following form of a geodesic equa- 
tion 



Then, a 'wave- function' is introduced, which is nothing but another 
expression of the complex action, and which reads 




(5) 



mc Vu = —duS. 



(6) 



^ = e' 



iS/So 



(8) 
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where is a constant, needed to ensure the correct dimensionahty 
of the equation and which can be written So ~ mcX. Therefore, the 
complex velocity field becomes 

v^ = aa^ini/>. (9) 

Introducing the 'wave-function' in the motion equation Eq. ((T]), one 
obtains 

(a''inv + ia") a„9,.inV' = o, (lo) 

which, after some algebra, takes the form 



, ,=0. (11) 

It can be therefore integrated as a Klein-Gordon equation, which 
reads 

d^d.xb + ^ = 0. (12) 

Such a construction can be generalized to the case when the full 
velocity field V + W is taken into account in the motion equation 
and in the covariant derivative. However, this leads to the same 
results as regards the final quantum mechanical- type equations [21] . 

Now, we come back to Eq. (0. The complex velocity is decom- 
posed into its real and imaginary parts, Vv = Vv— iUv, so that the 
geodesic equation becomes 

v -i(u'' ^^d''y^d^{v.-iU.) = o. (13) 

The real part of this equation takes the form of a relativistic Euler 
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equation which reads 

Equation (|9]) gives the imaginary part of the complex velocity field 
Vv, in terms of the wave function modulus, \fP, under the form 

= -\d,hi-/P. (15) 

Inserting this expression into the right-hand side of Eq. (|14p . one 
obtains, after some calculations, the following equation of dynamics 

m 



^ = — = - . (16) 

ds m 2 \ y^p ' 



Another way of writing Eq. p6|l is 



1 2, fd^d^VP 



v^d^v, = -\'d. y ^ j , (17) 

Now, the complex velocity, ~ V,^ — iU^, being linked to the 
'wave-function' by Eq. (|9]), its real part writes, in term of the phase. 



V. = -^d^S, (18) 



where S is the real part of the complex action 5. 

At the non-relativistic limit (c — > oo), the Dalembertian d'^d^ = 
d'^ /c^dt^ — A reduces to — A, the four-derivative d^, becomes —V and 
the fractal parameter A becomes 2V = caI^. Since in the quantum 



^Actually, the relativistic definition of A is: {d.^'^d^'^) = ^pAr^^^ds and the non- 
relativistic definition of "D is (d^±idi;±j) = ±2'DSijdt, where ryf" is the Minkowski metric 
tensor and Sij is the Kronecker symbol. Since r]^^ = —Sij and in proper time ds = cdt, 



14 



domain A — h/mc, we recover the non-relativistic standard 'quantum 
potential' Q = -h^ /2m A^/P/^, which can be interpreted as a 
byproduct of the underlying geometry of the fractal space [14lll5lfT9] . 
At this limit Eq. (|17p becomes 



and Eq. (|18|l becomes 



V-^, (20) 
m 



We recognize in Eq. (|19p the standard Euler equation with a 'quan- 
tum potential' obtained when performing a Madelung transforma- 
tion of the Schrodinger equation with no exterior potential and in 
Eq. (|20p the standard definition of the Madelung fluid velocity as 
a gradient of the phase of the wave-function. Note, however, that 
Eq. (|20|) is here derived, since Vi, has been defined from the very 
beginning as the real part of the geodesic mean velocity field. 

Finally, the continuity equation is recovered as the non-relativistic 
limit of the imaginary part of the Klein-Gordon equation. From the 
definition of the wave function given by Eq. (|8]), we obtain 

i; = VPe'^^^°, (21) 

where S is the real part of the complex action, S, and \/P comes 
from its pure imaginary part. Now, we calculate d^dfiiji which we 



at the non-relativistic limit A becomes cA = 2©. This can be readily verified by a 
dimensional analysis of the equations. 
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substitute in the Klein-Gordon equation (|12|) . Extracting the pure 
imaginary part and using Eq. (|18p . we obtain the continuity equation 
under the form 

a^(Pv;.) = o, (22) 

which, at the non-relativistic hmit, becomes 

^+V(PV)=0. (23) 

This imphes that the square of the wave-function modulus, P, can 
be regarded as linked to a genuine density, p, by mP — p, in the 
fluid interpretation of the space-time geodesies. Actually, in the 
scale relativity space-time approach, the probability density charac- 
terizes the potential geodesic bundle from the very beginning of the 
description. In this interpretation, the geodesic fluid is expected to 
be more concentrated somewhere and less concentrated elsewhere. 
It fills some regions and is nearly vanishing in others, as does a fluid. 
This behavior can be described by a presence probability density of 
the paths of the fluid elements and reduces to fluid motion equations 
in the classical domain where the classical part of the velocity fleld 
happens to dominate the fractal fluctuations [141 115j . 
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3 Quaternionic Klein-Gordon-type equa- 



tion 

In the following, we will consider the behavior of a system which 
is chaotic both in space and time. Such a property has often been 
ascribed to turbulence [22]. However, since a proper characterization 
of turbulence is much involved, we want to be cautious about the 
interpretation of our results which will be discussed in Sec. [5] 

First, we want to stress that, even if the mathematical construc- 
tion of the following is the same as the mathematical construction 
of the bi-quaternionic Klein-Gordon equation from which the Dirac 
equation is derived, its physical interpretation is different. Actually, 
we do not intend to claim that the macroscopic fluid whose features 
might emerge from our analysis will possess 'quantum' properties, 
neither would it be relativistic. This should be forbidden by the 
non-physical nature of the quaternionic Klein-Gordon-type equation 
which is unable to yield a Dirac equation for bi-spinors. Now, it is 
long known from experiment that the actual nature of the electron- 
positron couple is of the bi-spinor type. Therefore we consider the 
following quaternionic construction as typical of a representation of 
chaotic motion in space and time which is merely relevant at the 
non-relativistic limit and in the macroscopic domain. 

We want to stress once more that, in the scale relativity frame- 



17 



work, quantum-like equations have been proposed to act in the macro- 
scopic domain, providing therefore the constant A with a new status 
independent of the value of h. Actually, it is easy to verify that the 
whole mathematical structure of the theory does not rely on h and 
that it is preserved for any value of A which is introduced here as a 
transition length scale [TH]. This was first used for applications is- 
sued from the Schrodinger equation [5] , but it also obtains in the 
construction completed here. We therefore propose that Eqs. (|45p 
to H52p describe the behavior of a macroscopic non-relativistic fluid 
whose properties we examine below. 

Using the scale relativity method for deriving motion equations in 
a fractal space-time as recalled in Sec. [51 we introduce a quaternionic 
'wave function', which is defined as another expression for the 
quaternionic action and reads 



V-^a,V=-j9.S = -^V„ (24) 



where Vp is the differentiable part of the quaternionic velocity field 
and 5*0 /mc is a constant we rename A since it plays the same role of 
a transition scale as in the complex case. Ifere, ei is any of the three 
basis elements of the quaternionic algebra, {ei, 62, 63}, which satisfy 
61 = el = C3 = 616263 — —1. It is, for this algebra, an analog of 
the imaginary i for the complex number algebra. This is the reason 
why we use it to define the quaternionic 'wave function' from the 
quaternionic action. 
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We want also to stress that the 'wave function', -0, is here a 
mere mathematical tool used to express the action from which the 
velocity is defined. Since we do not consider this reasoning as valid 
to describe quantum mechanical processes we might have given to 
this tool any other different name. 

Since ■ip~^ij) — 1, which implies ip~^dfj,^p — — 9^V~^-V'i can 
write 



In the following, we will first compute the quaternionic compo- 
nents of the velocity field V,j under a parametric form to exhibit 
some of their physical properties, then we will employ the method of 
[19) to derive directly the fiuid motion equation from a free geodesic 
equation and finally we will extract the quaternionic imaginary parts 
of the Klein-Gordon-type equation to obtain conservation equations 
for the currents. 

For the sake of calculation simplification, we choose to write the 
quaternionic wave function as 



where g is any unimodular quaternion, and therefore such is g 
Inserting this expression in Eq. (|25p . we obtain 



= - ei A dfj,tp .tp. 



(25) 



^ = VPg-\ 



(26) 




(27) 
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Since g is unimodular, we can write it under the form 

g = exp(aei + /9e2 + 763), (28) 

where a, /3 and 7 are real fields. From this definition of g we obtain 
an expression for d^g.g~^ which we insert into Eq. (|27|l . After some 
calculations using the quaternionic algebra, we obtain a decomposi- 
tion of into its quaternionic components which are functions of 
the parameters a, /3 and 7. 

If we adopt the following notation for these components 

= -wom + ui^ei + W2pe2 + 113^63, (29) 

their parametric expressions are 

vo^ = X{d^a + sin 2/3 ^^7), (30) 

«i^ = ^9MlnP, (31) 

V2)i = A(sin2a9;j/3 + cos2Qcos2/?9p7), (32) 

U3p = — A(cos2a9p/3 — sin2Qcos2/3 9;j7). (33) 

It is easy, using this parametric representation, to check that the 
curls of the three fields vo,V2 and W3 are non zero and verify 

2 

dij.VOv — d^VOf^ = —{v2fiV3i, — V2vVifi,), (34) 

2 

dfiV2u - d„V2^ = -^(«3m"0i^ - «3i^"0m)i (35) 

2 

9fi"3i^ — d^V3f^ = —{vo^V2„ — Vo„V2^). (36) 
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As regards vi, which appears as a gradient of the logarithm of P, it 
wiU give rise below to the 'quantum' potential term in the motion 
equation. 

The quaternionic covariant total derivative applying here is given 
by Eq. ^ where the complex imaginary i is replaced by the quater- 
nionic imaginary ei as above (see [51 US] for an equivalent detailed 
derivation in the bi-quaternionic case). This gives 



Then, we insert this expression together with the decomposition of 
the velocity field V, given by Eq. p9p where vi has been replaced 
by its expression given by Eq. (jSTJ, into the free (inertial) motion 
equation ((Tjl. After some rearrangements, we obtain 



Now, for the sake of mathematical simplification, and as it has been 
done in Sec.[5J we derive the three other equations corresponding to 
the pure imaginary terms of the quaternionic motion equation from 
the Klein-Gordon-type equation, i. e., from the equation obtained 
when integrating Eq. ([7]). This yields already integrated equations 
under the form of continuity equations which are more easy to inter- 
pret. 

The quaternionic Klein-Gordon-type equation can be obtained 
with exactly the same reasoning employed to obtain it in the bi- 




(37) 
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quaternionic case [6l [13] . It reads 



(39) 



where i/' is the 'wave function' expressed by Eq. (|26p. 
We first expand 9'' 9^1/; as 

a^a^V = d^d^v^.g-^ + 2d''y^.d^g-^ + VPd''d^g-\ (40) 

Then, remarking that d^g~^ = ~'g~^{dy.g.g~^) and using Eqs. (|27p . 
(|29p and (|3ip . we obtain, after some algebra, an expression for d^dfiijj 
which we insert into Eq. (I39p where t/j has been replaced by its expres- 
sion given by Eq. (|26|) . Some final calculations give the Klein-Gordon 
equation under the form 



A'9^a^/P - Vp{vI + vI + vI) + Vp~ 2A d^'VPivof.ei (41) 
-"3^62 + ^'2^63) - \^{d^VQ^ei - d'^V3^e2 + a** 1)2^63) = 0. 



We extract each component of the quaternionic imaginary terms in 
the Ihs of Eq. (70) and equate them to zero. After some arrange- 
ments, we obtain three continuity equations which read 



(Pvo^) = 0, 



(42) 



d''{Pv2^) = 0, 



(43) 



a'^(p«3M) = 0. 



(44) 
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4 Non-relativistic limit 

At the non-relativistic limit, c — > oo, the time component of the 
quaternionic velocity is V" = c and the fractal parameter A be- 
comes 22? (see the footnote in Sec.[2ll. Moreover, since our purpose 
is to use the mathematical construction described here to apply it to 
the derivation of motion equations for fluids in the macroscopic do- 
main, this fractal parameter is no more linked to h but it is supposed 
to take macroscopic values. Therefore, Eq. (|38p . where we define the 
fluid density p = mP, becomes 

(I +^-^) --(I +--^) --(I +--^) - - -21.^^ (^) ' 

(45) 

which can be regarded as a motion equation for a 'three fluid' with 
velocity fleld V = {uo,i'2,W3}, where the ViS are vector fields in 
3-space, and with 'quantum' potential 

Q = 2rnD^^^. (46) 
VP 

Note that the sign of this potential is arbitrary since we could have 
chosen a reversed sign for the Ihs of Eq. (|45|l . We remark that the 
quantity T> has the dimension of a diffusion coefficient. Moreover, 
it is this parameter which determines the transition from chaotic to 
non-chaotic behavior. 

Now, we take the non-relativistic limit of Eqs. (|34|l to (|36|l and 
we obtain the following relations between the curls of each of the 
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velocity field component and the cross product of the other two: 

V X uo = ^ U2 X us, (47) 

V X U2 = ^ U3 X Uo, (48) 

V X ^3 = ^ t;o X i;2. (49) 
We recover here relations analogous to those found in [18) . with an 
additional intervention of the coefRcient T> which naturally arises in 
our formalism and gives the equations the correct dimension. They 
show that, in the most general configuration, the three fields have 
non zero curls which are each perpendicular to the other two fields 
and that these curl amplitudes are inversely proportional to T> and 
therefore to the strength of the potential. 

Now, we remark that the cases when Vi = Vj x u^, for i,j,k — 
0, 2, 3 and circular permutations of the indexes are particular so- 
lutions exhibiting the properties of Beltrami-type flows. In these 
cases, Eqs. (|47l) to (|49} give V x = 2?"^?;; for i = 0, 2, 3 and the 
ABC flow relations read Ai — BjCk — B^Cj, Bi — CjA^ — CkAj, 
Ci — AjBk — AkBj and six other equations with circular permu- 
tations of the indexes. We recognize the characteristic properties 
of Beltrami flows which are known to be chaotic [23]. Since such 
flows are particular cases of ours, this conflrms that our equations 
describe also some chaotic motion. However, Beltrami flows exhibit 
a vanishing velocity divergence which is not the case of our general 
solution. We can therefore conclude that our equations are able to 
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describe richer behaviors. 

Finally, we take the non-relativistic limit of Eqs. (|42p to (|44[) and 
obtain 

g+V(pi;o)=0, (50) 

^^+Vipv2)=0, (51) 

^+V{pv,)=0, (52) 

which are continuity equations for the three fluid 'currents'. Under 
this form, these equations apply to compressible fluids, with density 
p{t,x,y,z). For incompressible fluids, p = const., they become 

Vuo = 0, (53) 

V«2 = 0, (54) 
Vu3 = 0, (55) 
while the 'quantum' potential in Eq. (|45p vanishes and Eqs. (|47p to 



(|49p remain unchanged. 

Having found the general equations issued from our construction, 
we want now to discuss briefly some particular solutions which might 
be of interest in peculiar cases. We flrst equate the expressions of 
dp/dt extracted from Eqs. (|50p to (|52p and obtain 

V{pvo) = V{pv2) = V{pv3). (56) 



The first above equation gives 



V(«o - V2) ^ Vp 
Wo — f 2 p 
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(57) 



A particular solution of this equation can be obtained after a sepa- 
ration of the variables and reads 

voi - V2i = — , (58) 
P 

where i — x,y,z and the a^s are integration constants. Analogous 
solutions can be found for the couples {v2,V3} and {v3,vo}, which 
we write 

V2i — V-M = —, (59) 

p 

V3i -vo, = —. (60) 

p 

Eliminating p we obtain 

VOi — V2i V2i — V3i V3i — VQi , 

= = (61) 

ai bi Ci 

Since the a^'s, biS and c^'s are constants (and not fields), we remark 
that the ratios of the differences between the velocity fields are con- 
stant in every direction, which is a further constraint pertaining only 
to this particular solution. 

Let us stress, for completeness, that in the more particular case 
when ai — bi = d for all i, we can have two cases. Either vq = V2 
(equivalently vo = vs) and therefore Eq. (|45|l becomes 

where we have renamed v the field (equivalently 112). In this case 
the curl of v given by Eq. (|49|) (equivalently Eq. (|48|l ') vanishes and 
we recover a motion equation for an irrotational fiuid with velocity 
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V, subjected to the 'quantum potential' Q with a minus sign. This 
is analogous to Eq. p9p . i. e., the standard Madelung transformed 
Euler-type equation with no exterior potential, but with a reverse 
sign for the 'quantum potential'. 

The second case is obtained when V2 = V3. Then, Eq. (|45p be- 
comes 

(I H-.O.V) .0 - 2 (I + = -2V^V (^) . (63) 

Since «2 = V3, the curl of vo given by Eq. (|47|) vanishes. We are 
left with a two component velocity field V — {vo,V2} whose first 
component is a gradient and the second one is rotational. Note the 
appearance in Eq. (|63|) of a factor 2 in front of the second term in 
the Ihs which cannot be inserted in a redefinition of the velocity field 

Now we come back to the general equations obtained in this sec- 
tion. Remember that the chaotic to non-chaotic transition is char- 
acterized by the macroscopic parameter D. It was proposed, when 
applying the scale relativity tools to macroscopic systems, that D 
is a relative quantity specific of each system and that a general- 
ized Heisenberg relation could give an order of its magnitude as 
SxSv ~ 22?, where Sv is a velocity dispersion and Sx a position 
dispersion 24 . Now, it is well known, from experiment, that a tur- 
bulent behavior occurs in a fluid with kinematic viscosity coefficient 
V when its Reynolds number Re = VL/v, where V and L are char- 



27 



acteristic velocity and length scales of the flow field, becomes larger 
than some critical value Rec = VtLt/v. It would be therefore in- 
teresting to check experimentally if Re^ could be related to T> by 
21? = VtLt = RccV. If this condition appeared to be fulfilled, it 
might be hints that some kind of turbulent fiow in a fluid could 
be represented by a three component velocity field V = {vo,V2,vz} 
obeying a motion equation given by Eq. (I45|) and three continuity 
equations, (Eqs. (|50p to ((52])). The curls of the three vector fields, 
vo,V2,V3., would be non- vanishing and each of them would be parallel 
to the cross product of the other two fields. However, more investi- 
gations would be needed to confirm the turbulent behavior of a fiuid 
with such characteristics. 

It is here interesting to remark that the spontaneous appearance 
of vorticity at the macroscopic level for a chaotic behavior in space 
and time can be viewed as corresponding to the spontaneous ap- 
pearance of spin in quantum mechanics from the fractal feature of 
space-time jSj. Note however that both do not emerge at the same 
level of the construction since the spinor corresponds to the wave 
function ip and chaotic-vorticity to the derivative level expressed by 
Eq. ((251). 
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5 Discussion and conclusion 

We have chosen to study the chaotic behavior of fluid motion with 
the methods of scale relativity in a fractal space-time corresponding 
to a 'space-time chaos', since they have already proved to provide a 
good framework for the characterization of simpler fluid behaviors. 

Using the same line of reasoning as that developed in [19] . we 
have obtained a quaternionic Klein-Gordon-like equation, which, at 
the non-relativistic limit, yields a representation of a classical fluid 
behavior more complex than those previously obtained from a mere 
Schrodinger-type equation |15) . 

For the simplest case considered here, this solution is composed 
of a motion equation for a three component velocity field submitted 
to a 'quantum' potential, together with continuity equations for the 
three fluid 'currents'. The curls of each velocity component are non- 
vanishing and parallel to the cross product of the other two fields. We 
have also shown that particular solutions of our equations possess the 
physical properties of Beltrami fiows which are known as exhibiting 
chaotic behaviors. Since our general solution is more complex than 
such simple fiows, we propose it could be able to represent some kind 
of richer chaotic behavior, possibly turbulent. 

Now, in our model, the chaotic to non-chaotic transition is char- 
acterized by a macroscopic parameter T>, which has the dimension 
of a diffusion coefficient. If some kind of turbulent behavior hap- 
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pened to be reproduced by our solution, it might be linked natu- 
rally, in our framework, to some 'critical Reynolds number', Rcc, by 
21) = VtLt = Reef, where v is the fluid kinematic viscosity coefii- 
cient. 

We have thus all the tools to examine the way the transition 
from laminar to chaotic (possibly turbulent) behavior could occur 
in our representation. Actually, the extended chaotic behavior we 
propose can be viewed as the fundamental state emerging from the 
basic statement of our construction, i. e., the full relaxation of the 
differentiability assumption that leads to the breakings of the sym- 
metries ds <-> — ds and dx^ <-» — dx'^ . The chaotic time behavior 
is the degenerate state when only the symmetry <-» — dt is bro- 
ken. It yields, in the simplest incrtial case considered here, a linear 
macroscopic Schrodinger equation which gives an Euler equation for 
an irrotational fluid that is therefore unable to exhibit a chaotic be- 
havior. 

However, there is no a priori reason for the parameter "P; yielding 
the Schrodinger equation from which a laminar behavior occurs to 
have the same value as the parameter giving a Klein-Gordon 
equation from which a chaotic behavior can emerge, even for the 
same fluid. We might suspect that the length scale related to the 
former is of the order of the mean free path of the fluid molecules, 
while we know from experiment that the transition to chaos occurs 
at much larger scales. Therefore Vt should be considered as defining 
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the scales characterizing the transition from non-chaotic to chaotic 
behavior, and possibly, the critical Reynolds number. 

We have also shown that Vt appears both in the 'quantum' po- 
tential expression and in the curls of the velocity field components. 
This allows us to propose three means of testing our results. 

Since the strength of the 'quantum potential' is proportional to 
T)t , one might measure T>t in the chaotic phase of a given fluid sub- 
jected to such a potential. Then one could verify if the transition 
from a chaotic state to a laminar behavior occurs axound some 'crit- 
ical Reynolds number' corresponding to T>t. 

A second test proposal is linked to our result that the ampli- 
tude of the velocity component curls are inversely proportional to 
©. It is therefore a prediction of our representation that the greater 
the 'quantum' or dissipative potential applied to the fluid, the less 
whirling and therefore the less chaotic it might be. This seems consis- 
tent with the observation that, e. g., spontaneous turbulence disap- 
pears when viscous dissipation erases it. However, it would deserve 
a more exact experimental confirmation. 

A last test would be to examine if the velocity component curls 
are inversely proportional to some 'critical Reynolds number' multi- 
plied by the fluid kinematic coeflicient, which, as we have seen, might 
be a quantity proportional to T>t, in the somewhat turbulent case. 

Therefore we propose that our solution, which we currently present 
as a possible model for chaos (and maybe turbulence) since both its 
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theoretical foundations and its physical properties seem well adapted 
to such a fluid state, should be put to the tests above described before 
deciding if it can be accepted or rejected as such a representation. 
Moreover, these tests should provide a more thorough understand- 
ing of the properties exhibited by the model: e. g., the first test, 
by showing how the transition might occur around some 'critical 
Reynolds number'. 

Now, we want to stress that our results have been obtained from 
the simplest geodesic equation representative of inertial motion in 
the 'fractal space-time' with no additional force or potential nor vec- 
tor or other (e. g., tensor) field. We can therefore expect that, if this 
proposal was validated as a representation of turbulence, it should 
be designed to model some turbulence regime of the simplest kind. 
Now, depending on the equation of state of the fiuid, the 'quantum 
potential' can be regarded as a pressure potential or a mechanical 
stress tensor. Therefore, we suspect that more complex behaviors 
should emerge from a motion equation of the kind we have derived 
and to which new terms would be added that might be obtained from 
richer properties of the scale relativity construction (an example of 
such an enrichment is the passage from a simple Euler equation with 
a 'quantum' potential issued from a linear Schrodinger equation to 
a Navier-Stokes equation with pressure and viscosity terms corre- 
sponding to a non-linear Schrodinger equation with a vector field 
\15\). This will be the subject of future work. However, the three 
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fluid representation should still obtain in such enlarged cases since it 
derives from the quaternionic nature of the equations which is itself 
a product of the symmetry breakings. 

Note that our new representation is a generalization of previous 
proposals aimed at modeling vorticity and/or turbulence. 

Fluid vorticity first emerged in our context from a Madelung 
transformation of a Schrodinger equation where the complex wave 
function was subjected to a vectorial (electromagnetic-type) field 
[31 [TS] . However, the equations which were obtained involved only 
one rotational velocity field. 

A phenomenological two-fluid model of the truncated Euler equa- 
tions was proposed in [2S] as a turbulence model. Its results were 
shown to fit well the predictions of the Eddy-Damped Quasi-Normal 
Markovian (EDQNM) theory |26] which is known to well reproduce 
the dynamics of the truncated Euler equations [27] . Actually, the so- 
lutions of these truncated Euler equations have been shown to obey 
approximately K41 scaling [28) . However, the K41 scaling does not 
reproduce large amplitude events typical of intermittency which is 
one among the still misunderstood properties of turbulence [22] . 

Another well-known multi-fluid representation is the two fluid 
model for Helium II at non-zero temperature. It was developed by 
Landau [29] and involves a 'superfluid flow' and a 'normal flow' . 
However, it must be stressed that in such a model the fluid cannot be 
considered as an actual 'mixture' of two different liquids, in a classical 
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sense, since no friction nor momentum transfer occurs between the 
two parts. One can speak, in a certain sense, of the 'superfluid' and 
of the 'normal' parts of the fluid, but this does not mean that it can 
be actually separated into two such parts [30]. We suspect that this 
kind of feature might apply to our proposal where the 'true' nature 
of the three 'flows' might be a same kind of 'entanglement'. However, 
this also remains to be tested. 

Now we hope that, if it were validated by experiment, a three fluid 
representation emerging from the most basic and natural principles 
and equations of physics might somehow improve our knowledge of 
complex fluid motions. However, we are aware that solutions of the 
equations provided here might be hard to find, probably harder than 
for standard Euler and Navier-Stokes equations which are already ac- 
tual puzzles for mathematicians. Our hope is however that physical 
understanding might also gain from the study of the simplest special 
cases which should be easier to handle. 
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